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•  transmitted  region.  The  factor  DT  depends  mi  the  Incident  ware  front, 
tha  curvaturaa  of  tha  interface,  and  tha  relative  indlcaa  of  eha  two  aadla. 
Wa  give  explicit  matrix  foraulaa  for  calculating  DT,  illuatrata  ita  physical 
aignificance  via  examples,  and  point  out  an  arroneoua  aolutlon  in  tha 
raeant  lltaratura. 

Mast,  wa  atud^the  transmission  of  a  apharlcal  electromagnetic  wave 
through  a  dielectric  ahall.  Tha  two  aurfacaa  of  tha  ahall  ara  apharlcal 
(either  concave  or  convex),  and  thalr  cantara  are  arbitrarily  located  in 
relation  to  tha  aourca  point.  Tha  field  aolutlon  determined  by  tha 
geometrical  optica  theory  la  given  in  a  aimpla  cloaed  form.  Special 
attention  la  given  to  the  lane  affect  of  tha  dielectric  ahall  which 
converta  tha  incoming  apharlcal  pencil  into  a  fo cueing  pencil. 

Finally,  wa  present  axtanalva  numerical  results  on  arrays  of  point 
sources  covered  by  practical  radomas. 
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In  a  previous  report  (3],  details  of  the  analysis  of  a  curved  dielectric 
radem  using  geometrical  optica  ware  considered  for  a  point  source.  In  thla 
report,  ve  first  consider  the  transmission  of  a  spherical  or  plane  wave 
through  a  curved  dielectric  Interface.  The  transmitted  field  la  proportional 
to  the  product  of  the  conventional  Fresnel's  transmission  coefficient  and  a 
divergent*  factor  DF,  which  describes  the  cross-sectional  variation 
(convergence  or  divergence)  of  a  ray  pencil  aa  the  latter  propagates  In 
the  transmitted  region.  The  factor  DF  depends  on  the  incident  wavefront, 
the  curvatures  cf  the  interface,  and  the  relative  Indices  of  the  two  media. 

Ve  give  explicit  matrix  formulas  for  calculating  DF,  illustrate  its  physical 
significance  via  examples,  and  point  out  an  erroneous  nolutlrn  in  the 
recant  literature. 

Next,  we  study  the  transmission  of  a  spherical  electromagnetic  wave 
through  a  dielectric  shall.  The  two  surfaces  of  the  shall  are  spherical 
(either  concave  or  convex) ,  and  their  centers  are  arbitrarily  located  in 
relation  to  the  source  point.  The  field  solution  determined  by  the 
geometrical  optics  theory  is  given  in  a  simple  closed  form.  Special 
attention  la  given  to  the  Ians  affect  of  the  dielectric  shall  which 
converts  the  incoming  spherical  pencil  into  a  focusing  pencil. 

Final] v,  ve  present  extersive  numerical  results  on  arrays  of  point 
sources  covered  by  practical  radomea. 
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I .  INTRODUCTION 


Many  practical  antatmas  arc  covered  by  radomas,  whose  effects  on 
cha  ant  anna  radiation  ara  of  considerable  importance,  especially  in 
today's  hi^i. -performance  radar /communication  systems.  In  the  peat 
quarter  of  a  century,  savaral  standard  analyses  have  bean  devised  for 
analysing  redoes  affacta.  Nona  of  them  is  exact,  and  Improvements  ara 
alvaye  needed.  The  present  report  described  an  effort  in  this  direction. 

A  typical  radome  problem  may  be  stated  as  follows..  Let  an  aperture 
antenna  A,  for  Instance  a  horn,  a  slot,  or  a  conformal  array,  radiate  a 
known  field  ^(r)  in  free  space  (see  Figure  }).  A  protective  shield  or 
radome  I  la  pieced  around  antenna  A.  The  proolam  is  to  determine  the 
radiation  field  £  for  the  composite  structure,  i.e.,  the  antenna  A 
radiating  in  the  presence  of  the  radome.  This  problem  has  received  a  great 
deal  of  attention  from  many  researchers  during  the  last  two  decades,  and  a 
so-called  "best  available"  method  for  attacking  this  problem  appears  to 
have  emerged.  A  brief  description  of  this  method  is  given  b^low 

(a)  In  the  vicinity  of  I,  the  incident  field  I*  is  not  a  ray  field 
(locally  plane  wave).  To  cirausvent  this  difficulty,  let  ba  resolved 
into  a  spectrum  of  plane  waves,  namely, 


iVr) 


4  -4  i  W  •  r 

dk_  |  dky  W(k)  e1 


a.i) 


3(k)  -  j*  dx  j  dy  E*(r)  «~ik'r 

Here,  5  ■  (k  ,k  ,k  )  is  the  direction  of  propagation  of  the  plana-wave 
x  y  z 

spectral  component.  The  spectral  wave  number  in  the  z-direction, 


(1.2) 
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(a)  free  apace 


(b)  radome  environment 


Figure  1.  Antenna  A  and  radome  l. 


Figure  2.  Two  choices  of  incident  directions:  A1  and  PI 


k  •  -  (k^  +  k*) 

■  *  y 


(1.3) 


may  be  real  (homogenous  plane  vave)  or  imaginary  (inhomogenous  plane 
wave).  The  weighting  factor  $(k)  ia  the  amplitude  of  the  plane-wave 
spectral  component  propagating  in  the  direction  k. 

(b)  For  each  plane-wave  component,  a  tranomlsslon  coefficient  matrix 

Mi  ^ 

Tg(k)  for  a  flat  dielectric  slab  can  be  obtained  from  any  standard  text  on 
EM  theory.  The  subscript  taro  of  ?q  indicates  that  it  is  derived  from  the 
assumption  of  a  plane-wave  incident  field.  The  transmitted  field  %*(?)  at 
point  2  on  the  outer  surface  of  I  is  calculated  from  the  formula 


iC( 2 )  -  |  dkx  dky  T0(S>  ff(k)  elR‘r 


(1.4) 


(c)  Once  ?C(2)  1s  known  for  all  points  on  the  outer  surface  of  '£, 
equivalent  surfoce  current  sources  (1(2),  £(2))  can  be  determined.  The 
convolution  of  the  source  with  the  Green's  function  gives  the  desired 
radiation  field  which  is  expressible  as 


f 

t(r)  -  i  (G.  1  +  G,  it)  da 

!  11 


(1.5) 


outer  Z 


The  approach  described  above  is  of  course  theoretically  sound.  However, 
Its  faithful  execution  ia  impractical  because  of  the  extremely  laborious 
numerical  integrations  In  (1.4)  and  (1.5).  In  the  well-quoted  analyses 
by  Paris  [1]  and  Wu  and  Rudduck  [2],  the  numerical  integration  in  (1.4) 
is  avoided  by  approximating  the  transmitted  field  at  point  2  by 


ik  •  r 

Sc(2)  ;  f0(*Q)  3(£0)  «  0  r 


(1.6a) 


where  the  incident  direction  is  determined  by 


3 


kg  -  actual  ray  direction  Al,  or  the  direction  PI  of  the 

Poyntlng  vector  of  f1-  (Figure  2).  (1.6b) 

Note  that  the  approximation  in  (1.6)  is  to  deacribe  f*  by  a  plane  wave. 

Since  the  radoae  is  in  the  near  zona  of  the  antenna ,  this  plane  wave 
approximation  for  described  in  [1],  [2]  does  not  seen  to  be  a  good  one. 

In  the  present  report,  we  approach  the  redone  problem  from  a  different 
viewpoint.  Instead  of  deconposiug  the  incident  field  Into  a  plane  wave 
spectrum,  we  approximate  the  finite-sized  antenna  A  in  Figure  3a  by  an 
array  B  in  Figure  3b.  Each  element  in  array  B  radiates  a  spherical  wave. 
Those  spherical  wavs  constituents,  transmitting  through  the  radome  I,  are 
superimposed  to  give  rise  to  the  desired  radiatiou  field  f  in  the  far  zone. 
Thus,  the  key  step  in  the  present  approach  is  to  determine  the  transmission 
of  a  spherical  wave  through  a  curved  dielectric  shell. 

We  shall  apply  geometrical  optics  to  solve  the  transmission  through 
the  curved  radome.  Specifically,  we  first  study  in  detail  the  different 
aspects  of  refraction  at  a  single  curved  dielectric  interface;  second,  we 
analyze  the  behavior  of  a  spherical  dielectric  shell;  and  third,  we  apply 
our  theory  tc  various  practical  radomes. 
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Figure  3.  An  aperture  antenna  A  inside  a  redone  is  approximated 
by  an  array  B.  Each  point  source  in  array  B  radiates 
a  spherical  wave. 


Figure  4.  Transmission  through  a  dielectric  shell  due  to  Incidence 
from  a  point  source  at  P^. 


II.  DESCRIPTION  OF  PROBLEM 


The  geometry  of  th«  redom*  problem  under  consideration  la  sketched 
In  Figure  4.  A  point  source  at  Pq  produces  a  spherical  wave  which  goes 
through  a  curved  dielectric  shell  with  nonunifora  thickness.  Ray 
techniques  are  used  to  determine  the  field  at  point  on  a  given  surface 
outside  the  shell.  First,  let  us  describe  the  various  elements  Involved 
in  the  problem. 

Coordinate  Systems  and  Time  Convention.  The  main  coordinate  system 
is  the  rectangular  system  (x,y,z),  whose  origin  in  chosen  at  the  source 
point  Pq  and  the  z-coordinata  is  in  :'.e  direction  of  the  beam  maximum 
of  the  antenna.  Other  coordinate  systems  at  points  P^,  and  P^  along 
the  ray  are  defined  later.  The  field  is  time-harmonic  with  the  time 
factor  exp(+jut)  which  is  suppressed  throughout. 

Source .  We  assuma  that  the  source  has  a  well-defined  "phase  center" 

A  A  A 

at  point  Pq,  the  origin  of  the  coordinate  system  (x,y,z),  and  radiates 
a  spherical  wave  denoted  by  (i\  if*) .  If  the  antenna  is  an  array  of  point 
sources,  it  is  necessary  to  consider  each  elament  in  the  array  saparately 
and  superimpose  their  final  fields  et  the  observation  points. 

Dielectric  radons.  The  redone  is  e  dielectric  shell  with  nonuniform 
thickness  of  relative  dielectric  constant  •  e/c^  or  refraction  index 
n  *  */T^,  and  is  bounded  by  the  inner  end  outer  surfaces  and  Ij, 
respectively.  The  inner  surface  (near  the  source)  is  described  by 
the  equation: 

s  -  f^(x,y)  ,  for  e^  <  x  <  and  ^  <  v  <  d^  .  (2.1) 

The  outer  surface  i*  Rivea  by  the  equation: 
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t  -  f2(x,y) 


for  *2  <  *  *  bj  *°^  ej  <  7  * 


(2.2) 


> 


It  la  not  nacassary  to  know  tba  analytical  fora  of  tha  functions  f^x.y) 

anti  fjCx.y).  In  confutation,  only  a  aat  of  diacrata  data  points 

{x  ,y  ,f  }  with  n  ™  1,7,..., II  la  naadad  for  tha  daacription  of  f(f.  or  f,). 
n  n  n  *  * 

Thaaa  points  ara  flttad  by  a  cubic  spllns  vhlch  glvea  autoaatlcally  first 

and  second  partial  derivatives  of  f,  l.e.,  3f/3x,  ?f/3y,  32f/3x2,  32f/3x3y, 
2  2 

and  3  f/3y  .  Thera  are  two  requirements  for  tha  cubic-spline  fit: 

(1)  the  data  points  can  be  distributed  over  a  random  grid,  but  they  auat 
be  dense  enough  to  describe  thu  fine  details  of  or  Ej);  (11)  the 
domain  of  tha  data  points  (a  <  x  <  b  and  c  <  y  <  d)  must  be  somewhat 
greater  than  tha  area  of  C  In  which  tha  incident  ray  Is  expected  to 
intersect  the  radons . 

Observation  points.  Observation  point  P^  la  located  on  a  preapeclfled 
surface  C^,  which  can  be  either  one  of  the  following  two  types: 

(1)  Spherical  Z ^  with  center  at  PQ  and  an  infinitely  large  radius. 

In  this  case,  Pj  is  in  the  far  field,  and  the  field  at  ?3 
calculated  by  the  ray  technique  is  the  final  result. 

(11)  Planer  E,  which  is  just  outside  the  redome  and  normal  to  the 

(•axle.  In  this  case,  we  have  to  Integrate  the  field  on  l 3  to 
obtain  the  far  fiald. 

In  latar  calculations,  we  use  mostly  the  spherical  E^  in  (1). 
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III.  FIELD  SOLUTIOH  BY  GtOKmiCAL  OPTICS 


For  a  given  incident  field  (f^.t1)  generated  by  tha  source  at  point 
Fq  (Fig.  4),  tha  asymptotic  solution  of  tha  fiald  at  point  Pj  ia  deter¬ 
mined  using  geometrical  ope leu  [3),  [4].  The  method  of  solution  la 
daacrlbed  below. 

A.  Hat hod  of  So'vKloa 

Consider  a  ray  in  dlraetion  (9,4)  a* tending  from  tha  aourca  point  Pq 
ro  tha  point  P^  on  E^.  Tha  aourca  region  (Region  I)  is  homogeneous  and 
iaotroplc;  henca,  tw-  ray  la  a  straight  llua  along  tha  unit  vector  r^. 
First,  tha  distance  r^  la  found  and  tha  coordinates  of  point  P^  are 
detemined.  Than  tha  unit  vector  normal  to  tha  surface  E^  at  point  P, 

a  * 

ia  found  (Flgura  3) .  Tha  plana  of  vectors  r^  and  aatabliahaa  tha 
incident  plana.  Tha  angle  between  these  two  vectors  ia  tha  incident 
angle  ®*.  Using  Snail's  lav,  tha  refraction  angle  a*  is  obtained,  which 

A 

aatabliahaa  tha  direction  of  tha  transmitted  wave,  r^,  in  Region  II 

(dielectric) .  Tha  ray  in  Region  II  is  a  sti  tight  line  along  tha  unit 

*  ‘1*1*  *  *  * 
vector  r^7.  Three  coordinate  systems  (s^y^r^),  (u^,v^,N^),  and 

«t  *t  * 

fxl,yl,rl2)'  wlch  co*°°  origin  at  polut  P^,  tre  than  established. 

They  belong  to  tha  incident  ray,  tha  surface  E^,  and  tha  transmitted 
ray,  respectively. 

Tha  incident  fiald  is  split  into  a  normally  polarised  field 

(ij°,3jp),  (E  •vector  normal  to  tha  incident  plana  at  F^)  and  a  parallel 
polarised  field  ($*p,t*n).  Tha  transmitted  fiald  at  point  P^  is  obtained 
as  follows: 


•  r. 


1  l 


8 


IJP  -  2 


(3.1) 


SJQ  -  tp 


*  r 
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la  which  tB  *nd  tp  are  trunawlaalon  coefficients  for  the  no 
parallel  polarited  Helde.  respectively „ 


1  and 


n  2 

C1  •“ 


Y  coa  o 


t 


coa  a. 


al 


ci  “ 


pi 


nl 


Pi 


Yq  coa  0l 


coa 

_  1  CO*  °1 

2q  cob  a*  °  coa 


Z  coa  a. 


0.2) 


Y  «  —  •  /-2  «  — i- —  y  »  —  •  /—  •  nY  n  •  / —  •  /c~  . 

T0  ZQ  Jii0  120»  ‘  T  «  7w0  RT0  *  7*0  r 

Note  that  tha  aubacrlpt  1,  In  for  example,  algniflaa  tha  field 
evaluated  at  poiat  1. 

The  transmitted  field  at  la  Incident  on  I,  at  point  P^.  Coordlnatea 
of  this  point  cau  be  found  from  the  knowledge  of  the  coordlnatea  of  point 

A 

P,  and  the  transmitted  rav  direction  r,..  The  field  values,  In  going  fro* 

4  4  4 

P1  t0  V  undergo  some  change  which  la  dependant  on  tha  divergence  of  the 
ray  Thus ,  we  have 

-Jkr. 


•  <orl2). 


12 


(3.3) 


in  which  k  -  nk^  la  the  wave  ntaiber  la  the  dielectric  and  DF^2  la  th* 
divergence  factor  for  the  pencil  of  rays  travailing  fro*  to  P2  in  the 
dielectric.  It  la  given  in  (3],  aa 


nr  .  m  4.  II  a  II.  i-I/2 

DF12  *  ql  r12  (1  *  q2  r12) 


(3.4) 
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In  which  end  ara  the  principal  curvaturaa  (inverse  of  tha  radii 
of  curvatura)  for  tha  ray  pancil  in  Ration  XI.  Thay  art  found  fro*  tha 
curvatura  matrix  of  tha  tranamittad  ray  at  point  P^.  Tha  matrix  itaalf 
ia  found  fro*  a  foraula  involving  tha  curvatura  matrix  of  tha  incidant 
ray  and  that  of  tha  aurfaca  at  point  P^.  Tha  curvatura  matrix  of 
tha  tranamittad  ray  pancil  at  point  P^  la  alao  uaad  to  find  tha  curvatura 
aatrix  of  tha  ray  at  point  P2  incidant  upon  tha  aurfaca  Ij. 

Having  tha  flald  lncidaat  upon  I2  at  point  ?2,  tha  ray  dlractlon  r ^ . , 
and  its  curvatura  matrix,  uo  can  proeaad,  in  a  manna r  aimilar  to  tha 
transmisaion  through  £,,  to  find  tha  fl<ld  tranamittad  through  Ej  *t  P^ 
(Pigura  6).  Thua,  a  unit  vector  Hj,  normal  to  Ej  1*  obtained,  and 

A 

togathar  with  r^  dafinaa  tha  incidanca  plana  at  point  P^.  Tha  lncldanca 

t  X  A  A 

angle  (era  1b  than  calculated.  Again,  Snell'a  law  ia 

invoked  to  find  tha  refraction  angle  i*  at  Pj.  Thii  angle  ■pacifies  the 

ray  direction  in  Region  II.  (ou'.aid*  the  radoma).  Throe  coordinate 

“l  *1  *  ...  "t  't  * 

ayateaa  (x2,y2,r12),  and  (x^y-j.r^)  with  coaaon  origin  c: 

point  P2  are  then  introduced. 

Tha  field  ($2,ff*)  incident  upon  C2  at  Pj  ia  raaolved  into  parallel 
and  normally  polaritad  fialdc,  from  which  tha  tranamittad  fialda  are  found 
aa  follows: 


|tn  .n  |in 

K2  2  E2 


*? 


f0r23 


tp  Hln 

t2  Hj 


Pt? 


Z  Htn 
Z0M2 


'23 


(3.5) 
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in  which 


tn  - - L. 

*2  l  ♦  v 


n2 


n2 


1  ce*  a2 

n  1 

com  aj 


C2  ‘  — v 


p2 


p2 


cot  a. 


co*  a. 


0.6) 


Th*  (laid  at  obaarvatlon  point  P^  ia  than  found  fro*  th*  transmitted  fluid 


at  Pj,  auch  that, 


I,  -  (BF„)  . 


-JkOr23 


n 


(3.7a) 


in  vhich 


DP 


III 


-1/2, 


111 


23 


.-1/2 


(3.7b) 


aad  q*11  and  q^11  *ra  th*  principal  curvaturaa  of  th*  ray  pancil  in 
Region  III.  Thav  ara  obtained  fro*  tha  curvature  matrix  of  tha 
tranaaittad  ray  at  point  Pj.  This  matrix  ia  obtainad  fro*  a  formula 
alraady  mentioned  in  connactlon  with  trana*laaion  through  1^,  Tor 
a  typical  factor  in  Sqa.  (3.6)  and  (3.7b),  tha  following  square  root 
convention  la  uaad: 

♦|fl,  if  f  ia  real 

f  ■  1//1  +  qr  •  ■ 

+j|f|,  if  f  la  imaginary 


(3.8) 


It  should  ba  mentioned  hare  that  wa  have  ignored  omltlpla  reflections 
in  tha  dielectric  redone  throughout  our  analysis.  Details  of  th*  analysis 
stay  ba  f  ound  in  [  5 ) . 

We  now  summarise  th*  final  results  obtainad  so  f.«r.  The  point 
source  at  Pq  (Figure  4)  radiates  a  spherical  wave  described  by 


-iV  -  * 

.  1  —  [(?  coadcoed  -  Q  >104)1  «■  (F  aiode©#*  ♦  Q  co»4)y  -  F  •to4«] 
<f/ V  (3.9a) 


S1  «  t0r  *  I4 


(3.9b) 


vhorn  (r,8,4)  are  apherical  coordinatae  with  origin  >t  Fg.  Tha  pattern 

functions  F(0,4)  and  Q<8,4)  in  (3.9)  are  liven.  At  point  F2  (Fl*ur«  s), 

»1 

M  doooapoaa  the  field  Into  two  coap*nenta  1b  tha  dl raet Iona  of  (x^y^,  !.«• 

tj  -  P<v  VV  -  *  Voi  *  ‘1  '  (3  l0) 


At  cha  obaarvatlon  point  FJt  wa  axpreea  the  flald  aa  followa 

15  -  *  (il-y\)?2  .  h?  -  Va  *  *5  • 


(3.11) 


Xha  two  coaponanta  of  in  (3.11)  are  found  fro»  tha  aatrlx  ^nation 


r**6 

E3  *2 


*r\ 


-  (DF12)(OF23)a 


-Jk0(“l2*r23) 


tlt2(*2*xl)  ■tjtj^'y*)]  ^**1] 


-lcJt;(yJ««J)  cjt^yj-yj)  |«}-yj 


(3.12a) 


or  note  coapaccly, 


*  -3kO(“rI2  *  r23>«  1 

E.  •  (DF)  a  T  E2  . 


(3.12b) 
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to  (3.12),  a  la  cha  rafractloa  Indox  of  eha  dialactrlc,  la  eh* 
fraa-apaca  wav*  auabar  and  t",  t“,  t?  art  cha  ooraal  a ad  parallal 

traoMaiaaloa  coaffielants  at  points  and  Fj,  raapactivoly ,  at  givan 
In  (3.2)  and  (3. ft).  Tha  tvo  di/ar panes  factor*  ara  glvan  la  (3.4)  and 
(3.7b).  Thair  calculations  constitute  tha  aajor  effort  of  tha  preaant 
solution. 


IV.  REFRACTION  AT  A  CURVED  DIELECTRIC  INTERFACE 


A.  Introduction 

The  roftoction  st  a  dielectric  Interface  Ip  of  fundamental  importance 
in  electromagnetic  theory.  If  the  interface  is  arbitrarily  cur* td,  the 
only  available  solution  la  the  orr  derived  by  the  geometrical  optica  theory 
(CO).  Such  a  solution  conaiata  of  tvo  main  Ingredients:  the  wl 1-known 
Fresnel  formulae  for  the  tranamlaalon  and  reflection  coefficienta  (due  to 
A.  J.  Fresnel  in  1823)  and  a  ao-called  "divergence  factor  DF."  Surprisingly, 
the  aolo'lon  of  DF  vaa  derived  as  early  aa  1915  by  A.  Gull* trend  [6],  but 
its  application  vaa  nor.  widely  recognised  in  the  electromagnetic/optical 
coamunity  until  very  recently.  Ir.  1972,  Deachstpa  (31,  (4)  redr  rived 
GuUatrand'ii  result  by  naing  "curvature  matrices"  for  describing  curved 
aurfacea/wavefronta,  thus  resulting  in  greater  clarity  and  simpler 
rnmputaClona . 

In  this  chapter,  we  supplement  Deschampa'  results  by  giving 
explicit  formulas  for  calculating  various  curvature  matrices  and  by 
Illustrating  the  physical  significance  of  DF  via  analytical  and  numerical 
examples.  Another  motivation  for  the  present  work  la  to  compare  our 
solution  with  tha  one  described  by  SnyJer  and  Love  (71- [91  for  the  same 
problem.  It  is  shown  th«t  these  tvo  solutions  art?  not  in  agreement. 

We  believe  that  the  Snyder-Love  solution  is  incorrect. 
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B.  Final  Solution  for  the  Refracted  Fields 


We  begin  with  a  statement  of  the  problem.  Two  infinite  dielectric 
media  with  refraction  Indices  n.  and  n2  arc  separated  by  a  curved  inter¬ 
face  I  (Figure  7),  which  i9  described  by 

l:  z  ■  f (x,y)  .  (4.1) 


The  origin  of  the  (x,y,z)  coordinates  is  at  the  source  point  0  in  medium 
1.  The  source  emits  a  spherical  wave,  whose  electric  field  at  an  obser¬ 
vation  point  r  *  (r,0,4>)  is  given  by  [for  exp(Jut)  time  convention) 

~jklr 


e1*  -  i- 


[ep(e  ,<j>)  +  «Q(e,*)] 


(4.2) 


where  k^  *  2*/\^  -  n^aj/c,  and  (r,9,4i)  are  the  spherical  coordinates 

with  origin  at  0.  The  problem  at  hand  is  to  find  the  transmitted  field 

“*■£ 

E  at  a  typical  point  2  in  medium  2,  and  the  reflected  field  E  at  a 
typical  point  3  in  medium  1. 

We  attack  the  problem  by  the  geometrical  optics  theory  (GO)  [3],  [4j. 
Referring  to  Figure  7,  let  us  concentrate  on  a  typical  incident  ray  in  the 
direction  of  emanating  from  the  source  at  0.  The  "outward"  normal  to 
surface  Z  at  the  refraction  point  1  is  N.  The  plane  defined  by  the  ray 
01  and  N  is  the  piane  of  incidence.  With  respect  to  this  plane,  we  resolve 
the  incident  fi*>ld  E*  into  two  components:  perpendicular  component  E^ 


and  parallel  component  ft 

,i 


i 

II 


We  introduce  a  scalar  u  such  that 


,  for  perpendicular  polarization 


"ll 


for  parallel  polarization 


(4.3) 


•ft;  -ff 

Similar  decompositions  and  notations  apply  to  E  and  E  .  Then,  the  filial 


solution  derived  from  GO  has  the  followiug  form 


Figure  7.  Refraction  at  a  curved  dielectric  interface  Z 


X.V 


t  -jk2b  i 

u(2)  -  (DF)2  T  e  C  uX(l) 

-jk.c  . 

ur(3)  -  (DF)3  Re  1  u  (1)  . 


(4.4a) 


(4.4b) 


The  various  factors  in  (4.4)  are  explained  below:  T  and  R  are  the  well- 
known  Fresnel's  transmission  and  reflection  coefficients  (for  a  planar 
interface) ,  given  by 


1  +  Y  ' 


1  -  Y 
1  +  Y 


(4.5) 


where 


n(ccs  ci^/cos  a^)  ,  for  perpendicular  polarization 


jn  (cos  c^/cos  a^)  ,  for  parallel  polarization 


n  ■  (n2/n^)  ■  relative  refraction  index. 


The  incident  angle  and  transmitted  angle  are  related  by  the  Snell's  law 


sin  a_  ■  —  sin  cs, 
2  n  1 


(4.6) 


For  n  <  1,  a  critical  incident  angle  c~c  exists,  where 


sin  a  ■  n  ,  if  n  <  1 
c 


(4.7) 


If  <2^  >  ac,  <*2  defined  in  (4.6)  becomes  comply  end  the  simple  ray  picture 
shown  in  Figure  7  is  lost.  It  is  not  immediately  clear  how  the  present  CO 
solution  must  be  modified.  Therefore,  in  this  paper,  we  exclude  the  case 
°1  >  ac  w^en  n  <  The  factor  (DF,  in  (4.4a)  is  the  so-called  "divergence 
factor"  [3]  of  the  transmitted  ray  pencil  at  point  2  in  reference  to  point  1, 
It  is  given  by 


1 


1 


(4.8) 


(DFK  -  - -  —  . 

/I  +  (b/R21)"  /I  +  (b/R22) 

Here  *^22^  are  the  two  principal  radii  of  curvature  of  the  transmitted 

wavefront  passing  through  point  1.  The  sign  convention  of  (or  Rjj)  *-* 

as  follows:  is  positive  if  the  transmitted  rays  In  the  corresponding 

normal  section  are  divergent,  and  is  negative  if  the  transmitted  rays 

are  convergent.  The  square  roots  in  (4.8)  take  either  positive  real  or 

negative  imaginary  value.  Thus,  (DF^  is  positive  real  (no  focus  between 

points  1  and  2  on  the  transmitted  ray) ,  positive  .Imaginary  (one  focus 

between  1  and  2) ,  or  negative  real  (two  foci  between  1  and  2) .  The  factor 

(DF)^  in  (4.4b)  is  the  divergence  factor  of  the  reflected  ray  pencil  at 

point  3  in  reference  to  point  1.  It  is  given  by 

(DF)  -  - - - - - - - - - .  (4.9) 

A  +'  (c/R3iy  A  +  "(c/r32) 

The  determination  of  the  four  principal  radii  of  curvature  (R23 

is  the  key  to  the  present  problem.  In  Section  IV. C,  we  give  an  explicit, 

step-by-step  description  of  their  determination. 

In  summary,  for  the  refraction  problem  in  Figure  7,  the  final  solutions 
for  the  fields  of  the  transmitted  and  reflected  rays  are  given  in  (4.4). 

This  solution  is  based  on  GO.  It  is  valid  for  high  frequencies,  and  for 
all  cases,  except  when  total  reflection  occurs  (n  <  1  and  >  a^) , 
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C.  Calculation  of  Curvatures  of  Refracted  Wavefronts 

For  an  arbitrary  interface  l  and  an  arbitrary  Incident  ray  01 
(Figure  7),  the  calculation  of  the  four  radii  of  curvaturea  (R^ *^j2’R'31 **32^ 
is  not  a  staple  task.  In  thla  section,  we  present  a  systematic  and  explicit 
procedure  for  doing  this  calculation. 

Coordinate  systems  at  point  1.  Consider  a  ray  leaving  the  source  at  0 
in  the  direction  (9,4),  which  intersects  the  surface  I  described  in  (4.1) 
at  point  1.  The  distance  a  is  determined  froa  the  non-linear  equation 

a  cos  9  *  f(x  »  a  sin  6  cos  4,  y  »  a  sin  8  sin  4)  .  (4.10) 

The  unit  vector  in  the  direction  of  the  Incident  ray  is 

m  x  ain  9  cos  4  +  y  ain  0  ain  4  +  z  cob  9  .  (4.11) 

The  unit  normal  N  of  surface  Z  at  point  1  is 

n  -  j  (-fx  x  -  fy  y  +  h  (4.12) 

2  2  1/2 

where  A  M  +(1  +  f  +  £  )  ,  end  £  *  £or  example,  la  the  partial  derivative 

x  y  x 

of  f(x,v)  with  respect  to  x.  By  defining  4  positive,  we  have  chosen  N  in 

a  a 

the  forward  direction  with  respect  to  the  incident  rsy.  Vectors  x^  and  N 
define  the  plane  of  incidence.  At  point  3,  we  introduce  four  orthonormal 
base  vectors: 

MAN 

0tl,yl'Rl?  for  the  incident  rey  01 
(»2»Y2>z2^  ^or  th*  transmitted  rsy  12 

A  »  A 

(Xj,y3,*3)  for  the  reflected  rey  13 
(u,v,N)  for  the  surface  I 
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We  choose 


yl  "  y2  "  y3  *  V  “  *  *  *1 


(4.13) 


■  a  unit  vector  nerval  to  tha  plana  of  incldanca 


Than  it  follows 

u  *  v  *  N  ,  x."  y  *  *_  for  n  '  1,2,3  .  (4. J  4) 

(inn 

The  incident,  transmitted,  and  reflected  rays  are  respectively  in  the 
directions 

a  a  * 

u  ain  +  H  cos 

A  A  A 

*2  ■  u  sin  a2  +  N  cos  a2 

A  A  A 

-  u  sin  -  H  cos 

where 

sin  o2  •  n-1  sin  ,  0  <_  a, ,  a2  <_  w/2  .  (4.1Sd) 

Note  that,  because  of  the  particular  choice  in  (4.13),  both  end  are 
always  positive,  and  have  values  in  (0,»/2). 

Curvature  matrix  of  surface  £.  At  point  1  on  surface  E,  the  following 
two  vectors  lie  in  the  tangent  plane  of  the  surface: 


(4.15a) 

(4.15b) 

(4.15c) 


rlx  •  x  ♦  f x  *  (4.16a) 

*ly  "  y  +  fy  »  (4.16b) 
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where  (x,y,s)  are  evaluated  at  point  1.  With  reapect  to  thu  baae  vactora 
(r^.r^  ),  tha  curvature  matrix  of  E  la  given  by  [10] 


*e,G,  - 

f.F, 

f , E,  -  e.F, 

1  1 

1  1 

11  11 

«E  * 

1 

.2 

(4.17) 

4 

/lGl  - 

*1F1 

*1E1  “  flFl 

where 

A  - 

•►(1  ♦ 

2  7 

f  ♦  f 

!}l/2 

*  y 

1 

L  " 

I  ♦  f2 

-  f  f  ,  C, 

•1  + 

f2 

1 

X 

*  1 

x  y  ’  1 

y 

« 

i,  " 

-a”1  f 

f  -  -a"1  f 

1, 

a 

i 

o 

i 

>-* 

• 

1 

XX  ’ 

i  xy 

yy 

All  (x,y,x)'a 

are  evaluated  at  point  1.  Now  we 

tranafer 

tha  curvature 

matrix  with  reapect 

to  (r^.r^  )  to  that  ^th  reapect  to 

(u.v),  namely. 

o 

M 

f 

A'1 

Qz  A 

(4.18) 

where 

-► 

rlx 

•u 

r,  *v 
lx 

A  - 

• 

.V 

•u 

r,  *v 

iy  J 

It  nay  be  ahovn  [10]  that  a  principal  radlua  calculated  from  (4.17)  or  (4.18) 
haa  a  poaltlve  (negative)  algn  If  the  normal  aectlon  of  the  aurface  benda 
away  from  (toward)  tha  normal  N.  For  example.  If  E  la  a  aphere  with  radlua 
p  and  the  normal  N  polnta  away  from  the  aphere  center,  we  have 


(4.19) 
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asas.  *• 


-  ~ . 


Wa  not#  that  th«  prtitot  algn  cenmiloo  for  tho  aurfaca  curvatura  la  tha 
mm  aa  that  uaad  la  [4],  hut  oppoalta  to  that  la  (3),  (10). 

Curvatura  matrlcaa  of  wavafronta.  Tha  Incident  wavefront  paaalng 
through  point  1  la  apharlcal  with  a  radlua  a.  Thua,  Ita  curvatura  matrix 
Qx  with  raapact  to  baaa  vactora  or  any  othar  orthonormal  baaa 

vectors,  la 


(*.20) 


Tha  curvatura  matrlcaa  of  tha  tranamittad  and  raflactad  wavafronta  paaalng 
through  1  ara  axpraaaad  with  raapact  to  baaa  vactora  (x.,,y2),  and  (x^.y^), 
raapactivalv .  Thay  ara  danotad  by  Q2  and  Q^.  Tha  aolutlon  of  Q2  la  found 
from  tho  following  matrix  aquation  [4]: 


n  1,  B,  •  B*  ♦  (n  coa  «2  -  coa  a^)  QE  (4.21) 

whara 


*  A  A  A 

X  U  Xj^V 

coa  0 

* ! 

m 

L*i'u  yi*v 

0  1 

"a  A  A  a’ 

_ 

*2#v 

coa  a2  0 

B?  * 

m 

A  A  A  A 

y2*u  y2*v 

.  0  1 

The  aolutlon  of  la  found  from  tho  following  matrix  aquation 
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(4.22) 


8j  Qj  13  •  B*  Qx  Bt  -  2 (cos  Qt 

where 

-cos  0 

0  1 

Princlpsl  rsdll  of  curvature  of  rsfrsctsd  vs vs fronts .  Ones  matrices 
Q2  sod  Q,  srs  dsterained  frost  (4.21)  and  (4.22),  thsy  my  b«  dlagonallted 
In  s  standard  manner  to  find  thair  eigenvectors  (principal  directions  of 
tha  wavefront)  and  thair  eigenvalues  (principal  curvatures)  (10].  In 
particular,  the  principal  radii  of  the  transmitted  wavefront  (R2l‘R22^  *r* 
the  roots  of  the  following  quadratic  equation 

•^2  -  j[  (trace  Q2)  ♦  det.  Q2  •  0  .  (4.23) 

R 

If  Q2  in  (4.23)  is  replaced  by  Q^,  the  two  roots  are  the  radii  (R^R^) 
of  the  reflected  wavefront. 


*3  * 


Xj*u  Xj*v 


y3*u  y3*v 
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D.  Special  Cwi  Spherical  Inf  r fact 


To  illustrate  the  results  obt.alned  in  the  previous  cwo  sections, 
let  us  concentrate  on  a  special  case  In  which  the  interface  1  la  spherical 
with  radius  |p|,  as  shown  in  Figure  8.  Following  our  sign  convention,  the 
radius  of  curvature  of  I  is 

[♦|p |  ,  if  £  la  concave  when  viewed  fro*  tha  source  (Figure  8a) 


if  £  la  convex  (Figure  8b) 


(4.24) 


Without  loss  of  generality,  we  asaime  that  tha  incident  ray  fro*  the  source 
at  point  0  is  in  the  direction  (9 , 4  "  0).  The  plane  of  incidence  is  then 
the  x  -  s  plane.  Making  use  of  tha  foraulaa  in  Section  IV. C,  we  find  that 
the  principal  radii  of  the  transmitted  end  reflected  wavefronts  passing 
through  point  1  are 


R2I  "  co" 


a2>  [l 


2  ^  1  / 

cos  a,  ♦  —  (n  cos  a_  -  cos  3 
i  P  2 


i>] 


-1 
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1^1,  . 
—  ♦  —  (n  cos  -  cos  a, ) 
na  np  2  1 


-1 


•  I1  -  ' ]' 

|_a  p  cos  OjJ 

.  [a  .  2  c°*  *1} 

[a  o  J 


(4.25s) 

(4.25b) 

(4.26a) 

(4.26b) 


It  can  be  shown  that  Rj^  and  ere  the  radii  of  curvature  of  the  normal 
sections  in  the  x  -  s  plane  (plane  of  incidence) ,  whereas  R^j  end  R^2 
arc  those  in  the  orthogonal  directions.  Since  in  general  R^  1*  R^  end 
Rll  **  R32*  the  refracted  and  reflected  pencils  are  astigmatic. 
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(A)  Normal  lncldancii.  For  »  0 


(4.25)  and  (4  26)  become 

(4.27> 

(4.28) 


Thu* ,  for  normal  incidence,  both  r*fr*ct*d  pencils  h**.-a  spherical  wavefronts 
(no  longar  aatlgmatlc).  Th*  relation  in  (4.27)  may  b*  rearranged  to  te&d 


(4.29) 


which  Is  th*  wall-known  lana  equation  in  optica.  fSee  for  example 
Eq.  (40-14)*  p.  678  of  (11].)  Not*  tha  correspondin'  notations  used  in  [11] 
and  K*re:  n  •*  1,  n'  ■*  n,  a  -*  «,  a*  -*■  (-R^),  an<*  *  *  (-o).)  Th*  divergent 
iacid  -nt  pencil  from  a  point  sour>’o  la  converted  into  *  convergent  transmitted 
pencil  in  medium  2  when  R^  <  0.  This  occur*  whan 


« J 


>  0 


(4.30) 


If  0  >  0  (concave  dielectric  interface  shown  in  Figure  8a),  thie  la  possible 
if  n  «  (nj/n^)  <  i.  If  p  <  0  (convex  dielectric  Interface  shown  in  Figure  8b), 
this  is  possible  if  n  >  1. 

(8)  Comparison  with  Synder  and  Love’s  result.  In  a  recent  article  [7], 
Synder  and  Love  consider  the  problem  sketched  in  Figure  7  for  sn  incident 
plan*  wave  (source  distance  a  *  «  in  Figure  7).  Their  final  result  is  in 
disagreement  with  ours.  To  show  this  disagreement,  let  us  concentrate  or 
a  simple  caae  (Figure  9):  a  concave,  spherical,  dielectric  interface  la 
illuminated  by  au  incident  plane  wave  wt^ch  is  given  by 
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i  - 

k  (x,y,«)  -  t  a 


(4.31) 


The  problem  la  to  find  cha  high-frequency  solution  of  the  refracted  llelda 
along  the  y-axia.  Baaed  on  CO,  our  solution  la  given  In  (4,4),  (4.5),  (4.0), 
(4.9),  (4.27),  and  (4.28)  with  a  -*  *.  Using  the  coordinate  ayateai  in  Figure  9, 
the  tranamltred  and  reflected  flelda  are 

-Jk,y 


r(0.y,0)  -  t\- - - 


1  +  n 


L1  *  1  »  ! 

Er(0,y.0)  -  ■!_,  I  ;2|y|/o7J  Hr*  1 


v  >  0 


y  <  0 


(4.32a) 


(4, 32h) 


The  factors  in  [  ]  in  (4.32)  a’-e  divergence  factoiu.  The  intensity  or  power 
density  of  the  incident  field  on  the  central  ray  (along  x  *  t  •  0)  is  given  by 


I*  ■  Ke(y*(Ei  *  )}  -  (n^/120e)  watcs/m^ 


(4.33) 


which  is  independent  of  y.  This  is  due  to  rhe  fact  that  the  incident  field 
is  a  plane  wave  jnd  all  incident  rays  are  parallel.  The  intensity  of  the 
refracted  field  on  the  central  ray  does  vary  with  y,  namely. 


ic(v) 


_ 2_ 

1  +  nj 


l  ♦  [*-=-!]  f*j 

l  n  (Pj 


,  y  >  0 


(4. J4a) 


[1  -  n1 

2 

r  x 

U  +  nJ 

i  -  (2  |>  |/p). 

y  <  0 


(4.34b) 


At  the  focal  point  of  the  reflected  pencil  v  •  -p/2  in  medium  1,  the  intensity 
Ir  in  (4.14b)  predicted  by  the  present  GO  becomes  infinite  as  expected.  For 
the  same  problem  sketched  in  Figure  9,  Snyder  and  Love's  solution  is  given 
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in  Fqs.  f 29)  through  (34)  of  [7],  For  the  case  of  normal  incidence 
(c^  ■  0)  and  central  ray  (x  -  z  -  0) ,  Snyder  and  Love's  solution  reads 


?sL(0’y»0) " z  rhr  •  p  >s>  y  ” 0 

-  \  jk.,y 

V<J»y'0) " z  mre  •  y<0 


(4.35a) 

(4.35b) 


which  should  be  compared  with  our  solution  in  (4.32).  We  note  that 
(i)  divergence  factors  (DF)g  end  (DF)^  ate  missing  in  (4.35),  and  (ii)  the 
propagation  phase  factor  exp^-jk^y)  is  inissing  in  (4.35a).  Thus,  we  believe 
that  (4.35)  is  incorrect.  Furthermore,  for  each  Incident  ray  (fixed  a^), 
Snyder  and  Love  define  a  "power  transmission  coefficient  Tg,  "  by 
(Eq.  (35a)  of  [7]) 

f  \ 


SL 


(4.36) 


£ 

As  may  be  seen  from  (4.34h),  the  intensity  I  is,  in  general,  a  function 
* 

of  position  (x,y,z)t  because  of  the  divergence/convergence  of  the  reflected 
rav  pencil.  Then,  T^^  when  calculated  correctly  is  also  a  function  of 
position,  and  does  not  have  the  usual  significance  associated  with  the 
"power  transmission  coefficient." 
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•  **y.  -  • 


K.  Numeric*,!  Re&ultd  and  i)lscug9iou 

For  the  reflection  problem  sketched  in  Figure  1 ,  the  final  solution 
for  the  transaitted  field  uC  into  medium  2  is  giver  in  (A. 4a),  when 
the  incident  field  is  gi’»en  in  (4. 1).  In  this  section,  we  present  some 
numerical  results  for  uC  for  various  interfaces  and  source  locations. 

We  consider  three  types  of  interfaces:  the  spherical  incovfacc 
described  by 

(k/X1)  -  i  -  [1  -  (x2  +  v2)/xj]1/2  ,  ,4.37) 

the  paraboloidal  Interface  described  by 

(z/\)  •  (x2  +  y2)/2X2  ,  (4.38) 

end  the  hyperboloidal  interface  described  by 

(zAl)  “  I[l  +  2(x2  +  y2)/xJ]1/2  -  J  (4.39) 

where  X^  is  the  wavelength  in  medium  1  in  which  the  source  is  located , 

For  aasy  comparison,  we  have  chosen  the  above  Interfaces  such  that  they 
all  have  tne  mine  curvature  in  the  axial  direction  (Figure  1C).  There  are 
six  source  locations,  indicated  by  numerals  Inside  a  small  circle.  When 
the  source  is  at  location  1,  for  example,  medium  1  is  on  the  right  and 
medium  2  on  the  left,  and  the  interface  is  concave.  The  source  is  assumed 
to  be  y-polarized.  We  calculate  the  transmitted  field  in  the  E-plane 

A  A 

(plane  normal  to  x)  and  H-plane  (plane  normal  to  y) ,  In  these  two  planes, 
the  incident  field  Is  assumed  to  be 


2.0 


HYPERBOLOID  PARABOLOID 


Figure  10.  Three  dielectric  interfaces.  At  x  -  0,  ell  three  interfaces 
have  the  same  radius  of  curvature  of  1 


E-plane 


(4.40) 


e*(?) 


A 

e  1 

y  i 


H-plane 


Thus,  in  the  E-plane,  the  ^-vector  is  parallel  to  the  plane  of  Incidence; 
whereas  in  the  H-plane,  the  f-vector  is  perpendicular.  The  observation 
point  2  is  lu  medium  2  (Figure  7)  with  distance  b  •  (far  zone). 

We  calculate  the  normalized  far  field  defined  by 


ic(2) 

E- field  when  n^  i 

il(2) 

E-field  when  n^  ■  ^ 

Substitute  (4.4a)  and  (4.40)  into  (4.41).  Under  the  condition  b  ■+• 


(4.41) 


we 


have 


EN 


■"/vTi 


b  -*■  ® 


(4.42) 


Here  a  is  the  distance  between  the  source  and  the  Interface  along  the 
incident  ray,  and  T  is  the  Fresnel's  transmission  coefficient  given  in  (4.S). 
The  factor  ^R2lR22  tl'e  ra^^us  °f  the  Gaussian  curvature.  In  presenting 
the  numerical  results,  we  plot  EN  as  a  function  of  6,  where  0  is  the  polar 
angle  of  observation  point  2  measured  from  a  line  parallel  to  the  z-axls 
and  passing  through  the  source  point.  The  relative  index  n  ■  n2^ni 
is  always  set  at  2  (transmission  into  u  denser  me  Hum). 

Concave  spherical  Interface.  Figures  11  and  12  show  the  E-  and  H-plane 
.ar-field  pattern  EN  as  a  function  of  6.  Note  that  the  field  strengths 
Increase  as  the  source  moves  closer  to  the  Interface  (smaller  a) .  This 
is  mainly  due  to  the  fact  that  EN  is  inversely  i roportlonal  to  a,  according 
to  (4.41).  The  Gaussian  curvature  decreases  with  a„  but  not  enough 

to  offset  the  factor  (1/a)  in  (4.41).  For  source  3,  which  is  at  the  center 
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POLAR  ANGLE,  8 (DEG) 


of  the  spherical  lnurfiet,  oil  of  tha  incident  rays  art  no  rail  co  Cha 
Intar faca.  It  can  ba  shown  that  *21  •  “  *•  Thus,  EM  calculated 

fro*  (4.41)  is  equal  to  T,  which  is  0.667  for  tha  praaaat  caae  of  n  •  2. 

Of  particular  interest  is  the  H-plane  pattern  of  source  1  shown  in  Figure  12. 
Note  tha  Barked  asymmetry  in  tha  far-fleld  pattern  which  ia  due  to  tha 
asywsMtry  of  the  surface  with  respect  to  source  1.  Figure  13  shows  the 
variation  of  the  axial  far  field  whan  the  source  is  aoved  along  and 
parallel  to  tha  s-axla.  It  shows  clearly  the  Increase  of  the  field  as 
the  source  novas  closer  to  the  interface. 

Concave  paraboloidal  interface  (Figures  14  and  IS) .  The  far-fleld 
patterns  for  the  concave  paraboloid  are  quite  aieilar  to  those  for  the 
spherical  case.  However,  the  pattern  variations  are  sore  pronounced. 

Concave  hyperboloids!  interface  (Figures  16  end  17).  Note  that  the 
far-fleld  pattern  due  to  source  4  has  a  dip  Instead  of  a  peak  in  the 
axial  direction.  This  is  in  contrast  to  the  situations  in  Figures  11,  12, 

14,  and  15.  There  is  another  fact  worth  awntloning.  Because  of  the  choice 
of  the  same  axial  curvature  for  the  above  three  Interfaces,  the  axial  field 
is  the  sane  for  all.  interfaces  when  the  source  is  at  2,  3,  or  4.  However, 
for  source  1,  which  is  displaced  from  the  symmetry  axis,  the  normalised 
axial  field  EH(6  -  0)  increases  from  0.826  for  the  spherical  surface  to 
0.9S4  for  the  hyperboloid. 

Convex  interfaces.  The  E-  and  H- plane  far-field  patterns  for  a 
convex  sphere,  paraboloid  and  hyperboloid  are  shown  in  Figures  18  and  19 
for  source  locations  S  end  6.  The  source  locations  5  and  6  were  chosen 
based  on  (4.30).  Source  5  produces  e  divergent  axial  pencil  in  medium  2, 


: 
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-field  pattern  through  a  concave  paraboloid 


Figure  15.  Sew  a»  Figure  14,  except  for  B-plane  patten. 


bol oid 


POLAR  ANGLE,  0 (DEG' 

Figure  IV.  San*  an  Pigrre  it,  cuccpt  for  H-(  ine  pattern. 


09 


Figure  18.  E-plaiic  far-field  pattern  through  convex  interface*. 


Figure  19.  Sawe  is  Figure  12,  except  for  H-plMao  pattern 


whereas  source  6  produces  a  convergent  axial  pencil;  the  behavior  in  tliss. 
non-axial  direction  ia  governed  by  the  type  of  the  interface.  Thus,  .^s 
may  be  seen  from  figures  18  and  19,  the  far  field  In  the  axial  direction 
through  the  spherical  Interface  has  a  peak,  for  source  6  and  a  dip  for 
source  5.  This  is  also  the  case  for  the  paraboloid.  However,  this 
behavior  is  not  observed  in  the  hyperboloidal  pattern. 

For  all  the  convex  interfaces,  the  variation  of  EH  as  a  function  of 
6  in  (4.41)  is  predominantly  determined  by  the  radius  of  the  Gaussian 
curvature,  l^R21R22,  atu*  t0  a  ^easer  ex!:ent  by  T  or  a. 

Ray  Picture.  The  H-plane  pattern  due  to  source  6  for  a  convex 
sphere  is  given  in  Figure  19.  The  corresponding  ray  picture  is  shown 
in  Figure  20.  We  launch  6  rays  at  4°  apart  in  the  upper  half  x-z  plane 
(x  >  0).  The  transmitted  rays  are  first  convergent,  and  after  crossing 
the  caustic  surface,  become  divergent.  The  incident  rays  in  the  upper 
half  x-z  plane  within  a  20°  angle  give  rise  to  transmitted  rays  in  the 
lower  x-z  plane  (x  <  0)  within  a  13.5°  angle.  There  ^re  two  caustic 
surfaces  associated  with  the  transmitted  rays.  The  intersections  of 
the  caustic  surfaces  and  the  x-z  plane  are  indicated  by  crosses  and  dots. 
Similar  ray  pictures  can  be  drawn  for  the  other  cases  also. 
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F.  Conclusion 


lor  the  ref ruction  problem  sketched  in  Figure  1 ,  the  finsl  geometrlcal- 
optica  solutions  for  the  transmitted  field  and  the  reflected  field  are 
given  in  (4.4).  They  are  applicable  under  rather  general  conditions, 
namely,  tne  dielectric  interface  described  in  (4.1)  is  arbitrary,  and  the 
incident  field  in  (4.2)  from  a  point  source  is  arbitrary.  A  major  step 
in  calculating  these  3o  .utio.is  is  the  evaluation  of  the  divergence  factors 
in  (4.8)  aiH  (u.9),  which  involves  the  matrix  operation  described  by  (4.21) 
and  (4.22).  Strictly  apt iking,  the  present  solution  la  valid  in  the  high- 
frequency  limit  a  -*■  however,  practical  experience  has  shown  that 
solutions  of  the  present  type  ar<c  reasonably  accurate  as  long  as  the 
radii  of  curvutvre  of  the  dielectric  interface  are  in  the  order  of  a 
wavelength  or  i«ore. 
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V.  WAVE  TRANSMISSION  THROUGH  A  SPHERICAL  DIELECTRIC  SHELL 


A.  Introduction 

One  of  the  fundamental  problems  In  electromagnetic  theory  is  the 
transmission  of  a  spherical  wave  through  a  dielectric  shell.  This  problem 
has  numerous  applications  in  antenna  radomas,  electromagnetic  shielding,  and 
scattering.  It  appears  that  solutions  to  this  probleti  are  available  only 
for  the  special  case  where  the  shall  is  an  Infinite  dielectric  half  space. 
That  case  was  first  studied  by  Soamerfeld  in  1909,  whereas  later  research 
was  summarized  in  a  book  by  Brekhovskikh  (Chapter  IV  of  [12]).  In  this 
chapter,  we  consider  a  more  general  case,  namely,  ths  shell  has  two 
spherical  boundary  surfaces.  Unlike  the  Sommerf eld's  problem,  our  case 
does  include  the  effects  of  the  shell's  -.-.urvature  and  thickness.  Therefore, 
its  solution  should  be  of  more  practical  Interest. 

To  solve  our  problem  rigorously,  the  spherical  wave  expansion  may  be 
used.  However,  due  to  the  fact  that  the  source  location  and  the  two 
dielectric  surface  centers  do  not  coincide,  the  translatloual  addition 
theorem  for  vector  spherical  wave  function  [13]  must  be  used.  (Our  problem 
is  roughly  comparable  to  scattering  by  three  dielectric  spheres.)  This 
theorem  leads  to  a  complex  series,  which  makes  it  very  difficult  to  generate 
numerical  results.  In  this  paper,  we  use  the  geometrical  optics  theory  (Go) 
[4],  [14]  to  calculate  the  transmitted  field  in  the  problem  sketched  in 
Figure  21.  Such  a  solution,  though  only  approximately  valid  for  high 
frequencies,  i3  given  in  a  simple  closed  form.  Thus,  it  allows  us 
to  study  the  "cause  and  effect"  of  the  various  parameters  iu  a  convenient 
manner  and  gain  physical  insight. 


47 


B.  Solution  for  Concentric  Spherical  Shall 


Lot  us  consider  the  wavs  transmission  problem  in  Figure  21.  The 
boundary  surfacea  of  the  dielectric  shell  are  two  concentric  spheres 
with  radii  R^  and  R^  (both  positive),  and  with  a  common  center  at  Q. 

The  point  source  at  0  emits  a  spherical  wave  described  by  (for  exp  jut 
time  convention) 

-JV 

ll(r)  -  tP(e,d)e  ♦  Q(e,m)  .  (s.i) 

Here,  (r,9,9)  are  spherical  coordinates  of  r  with  origin  at  0.  Functions 
P  and  Q  describe  the  radiation  pattern  of  the  source.  The  wave  manber 
kQ  -  "  “<V0)1/2  **  thmt  *P*C**  >7  using  GO  [*],  [1*1, 

the  transmitted  field  at  an  observation  point  3  is  to  be  determined. 

Without  loss  of  generality,  we  assume  that  point  )  .  a  in  the  (x, a) -plane 
with  rectangular  coordinates  (x^,0,z^). 

A.  Ray  Tracing.  In  accordance  with  Snell's  lav,  we  trace  a  ray 
f ran  source  point  0  to  observation  point  3,  via  refraction  points  1  and 
2  (Figure  21).  Clearly,  all  four  points,  0  to  3,  lie  in  the  seme 
(x,c)-plane.  For  a  given  launching  angle  9  of  the  ray  and  the  distance 
c,  the  other  geometrical  parameters  can  be  determined  from  the  following 
relations: 

sin  a*  »  (d/R^Jain  9  ,  a  *  R^[sln(9  -  a*)]/sin  9  (5.2) 

sin  a*  •  n”*’  sin  a*  ,  sin  a^  ■  (Rj/R^)  sin  a^ 
k  *  ^2^*^°^®!  “  o|)]/sin  a*  ,  sin  •  n  sin 

x3  •  a  s:n  9  ♦  b  sin(9  -  aj  +  aj)  +  c  sin  (a*  +  a*  -  aj  -  a*  +  8) 
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•  cot  9  *•  b  cot  (•  -  4-  o£) 


♦  c  cot 


(a*  + 


o2  -  «1  - 


^2  *  «> 


Thut,  for  t  |lvtn  (t,c),  «•  c.to  determine  tht  position  of  point  3 
straightforwardly.  On  tht  othtr  hand,  when  point  3  it  given,  explicit 
formulae  do  not  txitt  for  determining  9.  Ons  hat  to  find  9  by  trial* 
and- error. 

B.  Field  on  tht  Bav.  Tht  prtttnt  vector  field  problem  can  be 
decomposed  into  tvo  tealar  ones:  one  with  the  electric  field  vector 
perpendicular  to  tht  plena  of  incidence  (f  ■  y  E^)  and  the  othtr  with 
the  electric  field  vector  parallel  to  the  plane  of  incidence  . 

The  final  tolutlon  for  the  trananitted  field  £*  at  point  3  derived  by 
GO  it  given  by 


ro> 


*i  «> 

ta  si  (i) 

-  (OF) 

05 

• 

T„  8^  (1) 

JkQ(nh+c) 


(3.3) 


The  varioua  factort  in  (S.3)  are  explained  belov:  Tj_  and  T^  are  the 
producta  of  che  transmission  coefficient!  at  points  1  and  2,  given  by 


r  ti 

-i 

r 

cot  a. 

i  co>  «■> 

.  .  *  - 

’  s 

n  i 

cot  a. 

cot  a. 

4 

(5.4a) 
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r  t*i 

.  cot  a. 

cot  a. 

•  *"  a 

n  i 

cot 

1  i 

cot  a2 

•»  « 

-1 


(5.4b) 


The  tvo  components  of  the  incident  field  are  calculated  fro*  (5.1): 


Ei-a)  -  ~ — 

o 

• 

♦ 

>-✓ 

cr 

-Jko» 

E||  (1)  -  1T— 

P(9,p  •  0)  . 

(5.5a) 

(5.5b) 
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The  derivation  of  the  divergence  factor  DF  In  (5.3)  requires  mm  effort. 


With  tha  details  given  la  Appendix  A,  tha  final  result  for  DF  reada 

DF  -  (1  +  •cub)*l/2  (1  ♦  *12b)‘l/2  (1  +  «21c)”1/2  (1  +  *22c)"1/2’  (5*6) 

At  point  1  or  2,  it  can  be  shown  that  the  principal  directions  of  tha 
transmitted  wavefront  are  precisely  the  two  directions  parallel  and 
perpendicular  to  the  plane  of  incidence.  Here  are  the  two 

principal  curvatures  of  the  transmitted  wavefront  at  point  1,  and 
^*21**22^  >r>  thos*  *c  point  2.  They  are  calculated  from  the  relatione 

<U  ■  (n  cos2  «i)-1jj  co.2  a*  +  i-(n  cos  a*  -  cos  a*)j  (5.7a) 


'l2  *  m  +  C<C0‘  “l  '  JT  “*  «1> 


(5.7b) 


(cos2  a*)”1  +  Kii)"1  n  co*2  °2  +  r“^co*  “  n  cos  1 


-lb-1 _ 2  i  A  1 


(5.7c) 


/v  .  -1.-1  1  .  t  i.. 

»ci2  •  n(b  v  ^ij)  +  jq(coe  «2  “  n  co*  ♦ 


(5 . 7d) 


The  sign  convention  of  *  is  aa  follows.  If  ic  is  positive  (negative), 

the  normal  section  of  the  wavefront  is  divergent  (convergent) .  For 

example,  if  the  transmitted  wavefront  at  point  1  is  the  same  as  the 

incident  spherical  wavefront,  we  have  <n  »  <12  «  +R~*\  For  a  typical 

factor  in  (5.6),  the  square  root  convention  is  • 

+|  f | ,  if  f  is  real 

f  -  (1  +  <b)'1/2  -  ■ 

<*j|f|,  if  f  is  Imaginary 

When  f  is  imaginary,  it  means  chat  the  ray  has  crossed  a  focus  01  the 


(5.8) 


ray  pencil.  The  (+j)  accounts  for  the  well-known  (v/2)  phase  retardation. 
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Several  general  co— ante  about  tha  solution  In  (5.3)  era  In  order: 

(1)  For  tha  case  where  a  total  reflection  occurs  (aj  or  a*  be canes  complex) , 
the  field  In  the  transaltted  region  Is  not  an  optical  field,  and  the 
present  ray  solution  (5.3)  Is  no  longer  valid.  (11)  Except  for  special 
cases,  e.g.,  normal  incidence  ■  0,  the  two  curvatures  (*21 ’*22^  °* 
the  transmitted  wavefront  Merging  from  the  dielectric  shell  are  not 
equal.  Thus,  the  transmitted  pencil  la  generally  astigmatic.  (Ill)  It 
Is  possible  that  and/or  <22  mr*  negative.  Then  the  divergent  pencil 
from  the  source  Is  transformed  Into  a  convergent  (focusing)  pencil  after 
propagating  through  the  dielectric  shell,  (iv)  The  solution  In  (5.3) 
remains  valid  for  more  geometries  than  the  one  shown  in  Figure  21.  This 
is  discussed  further  in  the  next  section. 

C.  Generalization  of  Final  Solution 

For  transmission  through  the  concentric  spherical  shell  in  Figure  21, 
the  final  solution  consists  of  two  parts: 

Part  A:  Ray  tracing  formulas  in  (5.2) 

Part  B:  Field  solution  in  (5.3)  through  (5.5) 

It  can  be  shown  that  Part  B  is  vslid  under  a  more  general  condition 
(Figure  22),  namely, 

(1)  The  centers  of  the  spherical  shells  Q1  and  Q2  need  not 

coincide,  as  long  as  the  four  points  (0,1, 2,3)  are  coplanar. 

(11)  The  surfaces  of  the  shell  can  be  either  concave  or  convex. 

Looking  from  the  source  side,  (or  R2)  is  positive  If  the 
surface  is  concave,  and  R^  Is  negative  if  the  surface  is 
convex . 
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(a)  R,  >0 ,  R2>0  (b)  R,  >0 ,  R*<0 


x-z  PLANE 


(c)  R,  <0 ,  Rt  >0  (d )  R,  <0 ,  Rt<0 


Figure  22.  General  conf igurations  where  the  field  solution  (2.3)  can  be 
used.  Points  (0.1. 2. 3)  must  be  coplartr. 


Condition  (i)  it  necessary  in  order  i or  the  acalar  transmit sion  coafflclanta 


and  V|j  in  (5.4)  to  ba  valid.  Whan  tha  four  points  (0,1, 2,3)  ara  not 
coplsnar,  tha  parpandicular  and  parr! lei  polarisations  are  nr  longer 
uncoupled.  Than  tha  scalar  or  T^  in  (5.3)  must  ba  roplacsd  by  a 
(2  x  2)  matrix.  Also,  (5.7)  bacomas  mors  complicated. 

For  tha  general  configurations  in  Figura  22,  tha  ray  tracing 
formulas  in  (5.2)  ara  not  valid.  However,  by  following  Snail’s  law, 
tha  ray  tracing,  avan  in  tha  moat  ganaral  situation,  is  concaptually  simple. 
Thus,  Instead  of  working  out  a  sat  of  ganaral  formulas,  wa  leave  it  to  tha 
individual  problems . 

D.  Axial  Incidence  on  Symmetrical  Shells 

To  study  tha  features  of  tha  present  ray  solution,  let  us  concentrate 
on  a  special  case,  whare  tha  four  points  (Qj.Qj.O.S)  ara  on  a  straight 
line  (Figura  23).  Than  tha  four  curvatures  in  (5.7)  reduce  to 


*  ■  k  .  L.  + 

11  12  na  n  Rx 


*22  "  *21 


nR^  ♦  an(n  -  1) 
bR^  ab(n  -  ?.)  na  R^ 


1  -  n 


(5.9a) 


- 2— ~  .  (5.9b) 

b  *  «-}  '2 


An  irterestlng  question  is  whan  does  or  <^2  become  nagstlva  (meaning 
a  convergent  pencil)?  This  is  answered  brlow: 

(1)  Negative  The  transmitted  pencil  inside  the  dielectric 

shell  is  a  convergent  one  whan  <  0  or 


R,  >  0 
1  -  n  1 


(5.10) 


If  R^  >  0  (concave  dielactric  interface  shown  in  Figure  24a),  this  is  possible 
if  n  <  1.  If  R^  <  0  (convex  dielectric  Interface  shown  in  Figure  24b), 


chit  la  possible  only  If  n  >  1.  Tha  dlatanca  bacvaeo  point  1  and  focal 
point  la 

(11)  Negative  «22.  Tha  transmitted  flald  In  tha  frae-space  region 
outaida  the  dielectric  ahell  la  convergent  whan  *c22  <  0.  Let  ua  concentrate 
on  a  special  caaa  In  which  the  thlcknaea  of  tha  dielectric  ahell  la  anaall 
so  that 

b  «  .  (5.11) 


Then  <22  in  (5.9b)  bee ones  approximately 


f  2  -  +  (n  -  1) 


1_  1_ 
Rl'  R2 


(5.12) 


which  la  the  well-known  equation  for  the  thin  Lena.  (See  for  example 
Eq.  (41-1),  p.  685  of  [11].  Note  the  corresponding  notation*  uaed  in  [ll] 
and  here:  a  •  -a,  a'  •  -*22,  Rj.  "  -Rj.  and  Rj  •  -Rj.)  The  condition 
for  a  negative  under  the  approximation  in  (5.11)  is 

*1*2 


(n  -  1)(R1  -  R2) 


>  0 


(5.13) 


(ill)  Far  Field:  If  the  observation  point  3  la  In  the  far  xone 


(<22c  >>  then  OF  in  (5.6)  becomes 
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U 

i  4  .  n  *b 

r2  *  h  i  ,  1 
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a  +  —  b  +  cj 
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^  p  p 

R1R2 
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OF 

The  first  factor  in  {  }  in  (5.14) is  the  divergence  factor  of  a  planar 
dielectric  slab  (R^.Rj  ■+  •).  Thus,  the  ratio  of  the  electric  field  at 
a  far-f laid  point  3  for  a  spherical  shell  and  chat  for  a  dielectric 
slob  is 


-1 

K5.14) 
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_ 1^(3)  for  spherical  shell! _ 

|EC(3)  for  a  slab  of  same  thickness 


1  +  (n  -  lj 


ab 

[R2  -  R, 

R-R» 

l  b 

-  1  +  - 
n 


fl  ^1 

|  T  * 

1  ’ 

«  J 

1  J 

-1 


(5.15) 


As  a  numerical  example,  consider  the  case  In  which  the  inaer  and  cuter 
dielectric  surfaces  are  concave  and  concentric  (Q^  ■  Qj)  with  (R^/b)  *  2. 

We  plot  n  as  a  function  of  (a/b)  for  n  ■  0.5  and  a  »  3  in  Figure  25.  We 
note  that  n  can  be  substantially  different  from  unity.  When  a  *  R.  ,  we 
have  n  ■  1.  Thus,  in  this  interesting  special  case,  the  axial  far  field 
through  a  concentric  spherical  shell  and  that  through  a  planar  slab 
become  the  same.  Another  interesting  special  case  occurs  when  DF  -*  ~. 

It  means  that  the  paraxial  rays  emerging  from  the  dielectric  shell 
(Figure  23)  are  parallel  to  the  axis  so  that  they  focus  at  the  far-field 
point  at  infinity.  From  (5.14),  it  is  shown  that  DF  “*»  if  k ^  ^  0  or 

.  R  -  R_  +  b(n  -  l)/n 

I  '  -  b(„  -  i)/„i  •  <5-I6> 

Under  the  thin- lens  approximation  b  '  0,  (5.16)  is  reduced  to  the  well- 
known  leneruker  *  s  equation  (see  F,q.  (41-2),  p.  685  of  [11]).  In  the 
antenna  radome  application,  (5.16)  is  useful  in  the  determination  of  the 
enhancement  of  the  antenna  main  beam. 

(iv)  Multiple  Refraction;  For  a  given  source  point  0  and  observation 
point  3  in  Figure  21  or  22,  we  can  trace  two  types  of  geometrical  optics  rays. 
The  first  type  is  the  direct  ray  from  0  to  j  without  going  through  internal 
reflections  in  the  dieleccric  shell.  Its  field  solution  is  given  in  (5.3) 
which,  of  course,  is  the  main  contribution.  The  second  type  contains 
rays  which  bounce  one  or  more  times  inside  the  shell  before  reaching  point  3. 
We  now  consider  the  contribution  of  such  multiply  refracted  rays.  For  the 
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axial  incidence  caae  (Figure  23) ,  let  ue  further  apccialize  tne  geoaetry 


T 


(5.20c) 


slab 


fn  ~  A2 

"  U  ♦  iJ 


•xp(-j 2kQnb) 


Comparing  (5.2Gb)  witf  (5.3.7),  wu  note  tha  efface  of  tha  multiply  rafractad 
raya  is  accounted  for  by  replacing  T  by  T  Whan  tha  condition  in 

(5.20a)  is  not  met,  va  must  evaluate  (5.19)  numerically.  Let  us  define 
an  error  term 


( 1^(3)  |  including  p  multiply  rafractad  raya 


*  100X.  (5.21) 


112*0)1  including  all  multiply  refracted  raya  J 
In  figure  27,  we  plot  5Q  (including  no  multiply  refracted  raya)  and  6^ 
(including  one  uuJtlply  isfracted  ray)  va.  u  for  b  ■  0.75  and 
(a  +  c)/b  -  A.  Several  observations  are  aade.  (a)  For  commonly  used 
values  of  n  (between  1  and  3),  the  error  ia  13X  or  leas  except  at 
resonances.  A  resonance  occurs  whan  all  multiply  refracted  rays  emerging 
from  the  slab  are  in  phase  with  the  primary  ray.  For  the  configuration 
ir>  Figure  26.  the  resonance  condition  la 


(nkQb/iT)  -  (2nb/XQ)  ■  a  positive  integer  .  (5.22) 

For  &n  obliquely  incident  ray  and/or  a  curved  slab,  the  condition  for 
resonance  is  rarely  satisfied.  Thus,  generally  speaking,  the  error  for 
neglecting  the  multiple  refraction  is  roughly  10X.  (b)  Errors  and  6^ 

have  about  the  same  order  of  magnitude.  Thus,  the  inclusion  of  the  first- 
order  multiply  refracted  ray  does  not  In  general  improve  the  accuracy  of 
the  solution. 
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E.  Nu—r teal  Results 


(a)  Radone.  Consider  a  dlalaccrlc  apharlcal  radome  with  an  inner 
radlua  20  A^,  a  unifora  thicknaaa  0,5  and  n  •  /5,  aa  shown  in  Figure  28. 
Tha  y-polarized  source  is  locatad  at  points  1,  3,  or  4,  and  its  radiation 


field  in  the  E-plane  is  given  by 

-JV 

^(r.e.d  ■  tr/2)  »  - — - —  [y  cos  (1.5  6)]  . 


(5.23) 


The  transaltted  field  is  calculated  fro*  (5.3).  For  an  observation  point 
in  the  E-plane  and  in  the  far-fleld  zone  (x  ■  0  plane  and  c  ■+  *»  in  Figure  21), 
we  may  rewrite  (5.3)  as 

"Jkor  . 

EC(r)  -  e__  ly  pC(Q)1  §  r  -  -  (5.24) 

where  r  is  the  distance  from  point  0  to  point  3  (Figure  21).  We  plot  pt(8) 
as  a  function  of  9  in  Figure  28.  Generally  speaking,  the  radoeie  Modifies 
the  radiation  field  gently,  as  expected. 

(b)  Lena  I  (Double  Concave).  Unlike  the  above  redone,  a  dielectric 
lens  nay  modify  the  Incident  field  drastically.  Let  ua  consider  Lens  I, 
drawn  approximately  to  scale  in  Figure  29.  The  source  is  2  XQ  away  from 
the  lens,  and  is  y-polarized.  Let  us  concentrate  on  the  field  in  the 
H-plane  (x-z  plane).  In  Figure  29,  we  launch  4  rays  2*  apart.  The  outside 
ray  (at  9  ■  6*)  suffers  total  reflection  at  the  second  face  of  the  lens, 
and  is  not  transmitted  into  the  free  space  region  (ve  Ignore  multiple 
refractions).  We  assume  the  incident  field  from  the  source  is  confined 
to  a  cone  (a  beam).  In  the  x-z  plane,  it  is  given  by 


if  9  <  9 1 
if  9  >  91 


(5.25) 
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0  ( DEG) 


1  dielectric  lens:  geometry 


In  tha  present  cat*  (Lana  I),  wa  chooaa  ei  -  5.5*.  After  transmission 
through  tha  lens,  tha  raya  become  more  divargant,  and  chay  are  apraad 
ovar  an  angular  raglon  of  about  0  <  53.5*,  which  la  conaidcrably  wider 
than  tha  incldant  angular  apraad  (6  <  5.5*).  At  an  obaarvatlon  point 
(r,0,$  ■  0),  wa  express  tha  transmitted  flald  aa 
■jk0r  . 

iC(r,6,0)  -  - — - —  (y  QC(r,0)]  .  (5.26) 

In  tha  abaanca  of  tha  lana,  QC  •  1  for  0  <  0*  and  QC  -  0  for  0  >  6*. 

t  3 

With  tha  Ians  praaent,  we  ploc  Q  ■' r , 0 )  aa  a  function  of  9  for  r  ■  2  *  10  A^ 

in  Figure  30,  Note  that  tha  transmitted  field  la  auch  weaker  (12%  or  less) 

than  the  incldant  field  beer  ise  of  tha  wider  apraad  of  tha  tranaaitted  rays. 

For  tha  present  case,  QC(r,6)  is  only  vary  weakly  dependant  on  r,  as 

long  aa  r  >  100  Xq.  Thus,  tha  tranaaitted  field  Cc  in  (5.26)  in  tha  far 

zone  is  approximately  a  spherical  wave  with  an  angular  pattern  Qt.  Every 

transmitted  ray  has  two  foci.  Their  distances  behind  the  second  face  E^ 

of  the  lens  are  (<0^)  *  and  (Kjj)  which  *ay  be  calculated  from  (5.7), 

In  particular,  (k^)  *  i*  ^ or  fhe  normal  section  of  the  wavefront  in  the 

plane  of  incidence  (x-z  plane),  whereas  (k^)”1  for  tf**c  in  th® 

perpendicular  plane  (defined  by  the  y-axis  and  the  ray  direction).  For 

the  third  ray  (incident  9  ■  4*)  in  Figure  29,  we  calculate  from  (5.7)  that 

(<21)”1  «  +0.21  XQ  ,  (k22)_1  -  +0.375  XQ.  (5.27) 

We  mark  the  position  of  the  focus  corresponding  to  by  a  cross  in 
Figure  29,  The  distance  between  A  and  the  cross  is  (tCj^,)”'.  Th*  trace 
of  the  two  sets  of  foci  is  shown  in  Figure  31.  They  are  curves  on  the 
two  caustic  surfaces  of  the  transmitted  wavefront  (intersection  of  caustic 
surfaces  and  the  x-z  plane). 
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Figure  10.  H-ntniu'  far  ti.it  forn  ilirimirli  lens  1  and  lens  II 


(c>  Lens  II  (Doubls  convtx) .  Thv  geometry  of  Lena  II  an d  Its  rsy 
picture  sr«  shown  in  Figure  32.  Ths  incident  field  is  given  in  (5.25)  with 
6*  “  10*.  Ths  transmitted  field  in  rhe  H-plene  is  expressed  in  (5.26),  where 
QC(r,6)  is  age in  very  weekly  dependent  on  r  in  the  far  tone  and  is  plotted 
in  Figure  30.  We  note  that  the  incident  10*-beam  la  now  spread  into  a 
67s* bean  after  traneaisaion  through  Lens  II.  The  peak  value  of  the 
transmitted  field  is  about  24X  of  the  incident  field  at  the  same  far-field 
location.  The  caustic  curves  are  shown  in  Figure  31. 

(d)  Lens  III  (Convexo-Concave) .  The  geometry  of  Lena  III  is  shown 
in  Figure  33.  The  radii  (R^.R^)  of  the  lens  and  the  source  distance 
satisfy  the  lensmaker'a  equation  in  (5.16),  so  that  the  transmitted  rays 
near  the  axial  direction  are  almost  parallel  and  focus  at  a  point  at 
infinity  in  the  axial  direction.  The  Incident  field  is  given  in  (5.25) 
with  e1  »  12*.  The  H-plane  transmitted  field  is  expressed  in  (5.26), 
where  Qt(r,9)  varies  drastically  from  the  near  field  rone  up  to  r  '  100  Xq 
as  seen  from  Figure  34.  Beyond  r  '  100  X^,  the  beam  becomes  narrower,  and 
the  peak  becomes  higher  with  the  increase  of  r.  It  is  well-known  that 
the  exact  value  of  the  beam's  peak  (on  a  caustic  surface)  cannot  be 
predicted  by  the  present  geometrical  optics  theory.  It  can  be  calculated 
from,  for  e.g.,  the  Huygens-Green  formula  described  in  p.  107  of  [15]. 

The  caustic  curves  are  shown  in  Figure  35. 
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Figure  32.  Double-convex  spherical  dielectric  lens:  geo^try  sod  ray  picture. 


0  (DEG) 

Figure  ?4.  H-plane  far  field  pattern  through  lens  III 


F.  Conclusions 


For  a  given  incident  field  in  (5.1),  a  geometrical  optics  solution  of 
the  transmitted  field  through  a  spherical  dielectric  shell  (Figure  22)  is 
given  in  (5.3)  through  (5.7).  The  present  solution  is  an  approximated  one. 
We  svomnarize  its  validity  and  limitations  below 

(i)  Because  ray  techniques  are  being  employed  for  the  analysis 
presented  herein,  the  radii  of  curvature  should  be  large 
in  terms  of  the  incident  wavelength  in  order  for  the  results 
to  be  valid.  As  in  all  ray-optical  solutions,  our  solution 
fails  when  the  observation  point  is  near  the  caustic  surfaces 
of  the  transmitted  wavefront. 

(ii)  The  solution  given  in  (5.3)  and  (5.7)  is  valid  only  for 
spherical  shells,  and  when  four  points  (0,1, 2, 3,)  in 
Figure  22  are  coplanar.  These  two  limitations  can  be  relaxed. 

In  fact,  an  explicit  GO  solution,  (which  is  very  similar  to 
the  present  one,)  for  an  arbitrarily  curbed  shell,  and 
arbitrary  source  and  observation  points  is  given  in  [5]. 

(iii)  Cur  solution  includes  only  the  contribution  from  the  first- 
ordev  refracted  ray  (a  direct  ray  from  the  source  to  the 
observation  point  withcvt  internal  refractions  inside  the 
dielectric  shell).  It  has  been  estimated  that,  except  at 
"resonances"  which  rarely  occur  for  curved  shells,  the 
error  of  neglecting  higher-order  refracted  rays  is  roughly 
102  for  n  »  /c^  in  the  range  1  to  3.  It  should  be  pointed 
out  that  the  higher-order  refracted  rays  can  be  calculated 
in  exactly  the  same  manner  as  the  first-order  refracted 
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ray  by  using  the  formulas  derived  in  this  work  and  in  [5]. 

It  is  only  a  matter  of  bookkeeping  and  computer  time  in 
doing  the  numerical  calculations. 

(iv)  For  transmission  through  a  curved  dielectric  shell,  the 

major  field  contribution,  of  course,  comes  from  the  refracted 
rays  which  are  studied  in  this  paper.  However  other  diffraction 
processes  exist  which  may  not  be  conveniently  fitted  into  ray 
descriptions,  and  their  contribution  may  be  significant  under 
certain  conditions.  An  example  is  the  modal  fields  guided  by 
the  curved  dielectric  shell.  Efforts  to  understand  its  behavior 
in  canonical  problems  have  been  initiated  [16]. 


VI.  ARRAY  COVERED  BY  RADOME 


In  this  chapter,  we  study  the  radiation  pattern  of  a  linear  array  of 
eleven  point  sources  enclosed  In  a  redone.  The  array  elements  are  arranged 
symnetrlcally  along  the  x-axia  with  a  half -wave length  spacing  as  shown  in 


Figure  36.  The  element  pattern  is  assumed  to  be  of  the  form 


I1  « 


-iV 


[sin  4  9  +  cos  ^  >]  cos  6 


(6.1) 


To  calculate  the  array  pattern,  rays  are  traced  from  each  source  through 
the  radome  as  discussed  in  Chapter  I.  When  the  ray  tracing  is  done,  there 
are  two  options  available  in  calculating  the  far  field:  (a)  we  can  calculate 
the  transmitted  field  E^  over  a  planar  surface  just  outside  the  radome 
and  then  integrate  this  field  to  get  the  far-field  pattern;  (b)  we  can  push 
the  surface  to  infinity  and  calculate  £!j  on  this  surface  directly,  thus 
avoiding  the  integration  step-this  Is  known  as  the  direct  ray  method. 

Though,  in  general,  both  methods  do  give  identical  results,  the  second 
method  becomes  less  accurate  if  the  far-field  point  happens  to  be  in  the 
vicinity  of  the  caustic  surfaces  of  the  transmitted  wavefront.  However, 
one  does  not  usually  confront  such  a  situation  in  practice.  The  details 
of  the  calculations  using  these  two  methods  may  be  found  in  [5]. 

Once  the  far-field  pattern  of  each  source  is  obtained,  the  sum  pattern 
is  obtained  by  adding  the  patterns  of  the  individual  elements.  A  progressive 
phase  shift  is  introduced  among  the  elements  to  enable  beam  scanning.  In 
the  results  to  follow,  we  present  the  H-plane  (x-z  plane)  sum  patterns.  Thus, 
restricting  ourselves  to  this  plane,  the  incident  field  in  (6.1)  reduces  to 
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ccs  8 


(6.2) 


I1  -i 


The  total  transmitted  field,  summed  over  all  the  sources,  may  then  be 
represented  by 


(6)  .  e _ y6  igt  i  #j (n-1)  (kdsine+a) 

sum  r  (8)  n 

n-4 


(6.3) 


where  |E  (e)|n  represents  the  magnitude  and  9^(6),  the  phase,  of  the  trans¬ 
mitted  field  due  to  the  nC^  source,  d  in  the  element  spacing,  and  a  Is 
the  inter-element  phase  shift. 

Array  patterns  were  calculated  for  four  types  of  radomes  as  specified 
in  Table  I.  Two  of  the  radomes  were  paraboloids  and  two  were  tangent 
ogives . 

(a)  Radoae  A:  This  is  a  paraboloidal  radome  cf  relative  dielectric 
constant  2.3.  The  inner  and  outer  surfaces  of  this  radome  were  generated 
by 

zAq  -  50  -  (x2  +  y2)  /  (8Xq)  (6.4) 

and 


z/\Q  -  50.25 


(x2  +  y2)/(8.16X2)  , 


(6.5) 


respectively.  The  geometry  of  this  radome,  along  with  the  eleven  element 
array,  are  shown  in  Figure  36.  The  H-plsne  element  patterns  through  the 
radome  are  shown  in  Figure  37.  The  H-plane  sum  patterns  are  shown  in 
Figures  58-42  for  different  scan  angles  (the  scan  angles  Indicated  in  the 
figures  are  the  values  for  the  beam  without  the  radome,  measured  from 
the  z-axls). 
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Figure  AO,  Seme  *s  Figure  38,  except  for  «cen  angle 


POLAR  ANGLE  ,  9  (DEGREES) 

Figure  41.  Seat  *•  Figure  38,  exc.«pc  for  scan  angle  »  -19.77*. 
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POLAR  ANGLE  ,  0  (DEGREES) 


Figure  42.  Same  ae  Figure  38,  excepc  for  scan  angle 


SAHPE 

e 

r 

LENGTH 

(External) 

DIAMETER 

(External) 

Paraboloid 

2  5 

50.25  xj 

40.50  XQ 

50.50  XQ 

40  XQ 

Tangent 

5.7 

14.6  :.Q 

6.7  XQ 

Ogive 


WALL 

THICKNESS 


_ 

y2-o 

0.233  XQ 


Figure  43  a how*  the  redone  bores ight  error  for  this  radome .  The  bora- 
sight  error  is  the  dlfrereoce  between  the  beam  maxima  with  end  without 
the  radome  ( (beam  aaxiaua  with  redome)  -  (bees  mexiauu  without  redome)  ] . 

It  may  alao  be  observed  from  Figures  38-42  that,  as  the  scan  angle 
increases ,  the  gain  decreases  because  of  the  decrease  in  the  effective 
aperture  with  increasing  scan  angle.  Figure  44  shows  this  gain  variation. 
Observe  that  in  the  presence  of  the  radome  the  maximum  field  strength  drops 
ouch  faster.  Also,  around  the  tip  region  where  the  curvature  is  the  highest, 
the  maximum  field  strength  decreases  faster,  and  then  levels  off  at  higher 
scan  angles. 

(b)  Radome  B:  This  is  also  a  paraboloidal  radome  of  -  2.S,  but 

the  outer  surface  in  this  radome  was  configured  so  as  to  make  the  curvature 
in  the  tip  region  very  close  to  the  critical  curvature  (see  Figure  30,  [3]), 
thus  resulting  in  a  greater  focuselng  action  in  the  x-dlractlon.  The  inner 
surface  was  obtained  from  (6.4)  and  the  outer  surface  was  obtained  from 

r/XQ  -  50. 5G  -  (x2  +  y2>/(7.925*J>  •  (6.6) 

The  geometry  of  the  radoce  is  shown  in  Figure  45.  Observe  that  the  thickness 
is  maximum  along  the  t-axis  and  reduces  gradually  to  xero  thickness  at  the 

base. 

The  H-plane  sum  patterns  are  shown  in  Figures  46-49,  and  the  bores ight 
error  curve  is  presented  in  Figure  50.  The  focussing  action  of  this 
radome  may  be  clearly  observed  in  Figure  46  -  the  axial  field  strength  is 
1.4  tines  as  large  as  that  without  the  radome. 


SCAN  ANGLE  ( DEGREES ) 


Figure  43.  Soreotghe  error  for  radow*  A. 


•i 


POLAR  ANGLE  ,  8  (DEGREES) 


Figure  46.  H-.  lane  sue  pattern  through  radoae  B.  Beaa  at  0*. 


90RESIGHT  ERROR  ( OEGREES ) 


(c)  Hadone  C  wd  Pi  The  a  a  art  tangent  ogive  rtdoMi  used  in  niaallaa. 
Tha  gecaatrlcal  configuration  la  givan  in  Figure  SI.  Radon  C  ha  a  -  S.7 

stui  radon  D  c  •  9.3.  It  aay  ba  observed  from  Flgura  51  that  the  radon 
ha a  a  discontinuity  In  tha  a-dlrtction.  In  practice,  thla  discontinuous  tip 
ration  io  raplacad  by  a  aaooth  curved  aurfaca.  However,  in  our  calculation, 
cMi  ■edification  waa  not  lntroducad.  Tha  (laid  along  tha  t-axla  waa 
calculated  by  Interpolation. 

Tha  K-pJana  sua  patterns  for  radon  C  ara  praaantad  in  Figures  52-57, 
anti  tha  corresponding  boraalght  error  curve  la  shown  in  Flgura  Tha 

pattern  fot  radoui  D  ara  in  Figurca  59-64.  The  Soraaight  error  curve  la 
ah  own  in  Flgura  65.  It  nay  ba  obaarvad  tram  thaaa  figuraa  that  tha  pattern 
for  anil  acan  anglaa  ara  not  generally  good.  Ala©,  tha  boraalght  error 
at  anil  acan  anglaa  is  quite  high.  Thta  behavior  la  not  uaually  obaarvad 
In  practical  radon a  for  two  reasons; 

(1)  In  practice,  the  entamu.  «  uaually  t  reflector  ontenna  or  a 
planar  array  of  considerably  large  nuatber  of  elannts,  vharaaa 
In  our  nodal,  we  have  only  elnvtu  point  aourcas  VQ/2  apart. 

If  ona  includes  aufflclant  mnaboi  of  point  eourcea  ao  aa  to 
faithfully  represent  a  givan  practical  antenna,  there  will  ba 
conaldertble  changes  In  the  eua  patterns  and  tha  boraalght  arror 
curves. 

(11)  Tha  aacond  reason  has  aooathlug  to  do  vlth  tha  tip  raglon. 

Flgura  64  shown  a  rengent  ogive  and  a  paraboloidal  radon  of 
the  ease  alee.  It  ny  ba  observed  that  tha  tip  raglon  in  tha 
rangact  ogive  le  very  eteep  whereat  in  the  paraboloid  tha 
surface  changes  acre  onoothly. 

In  ganeral,  the  patterns  lap rove  at  larger  scan  angles,  aa  dost  tha 


boraalght  error. 
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turn 


POLAR  ANGLE,*  (DEGREES) 

Figure  57.  Sim  as  Figure  52,  axe apt  for  scan  angle  -  -27.64*. 


103 


i aa&aflaL3: 


BORESIGHT  ERROR  ( DEGREES ) 


-30  -20  -10 


0 


10  20  30 


POLAR  ANGLE  ,  8  (DEGREES) 

Flgura  59.  H-plana  pat tarn  through  rado*a  D.  Baaa  at  0*. 
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SCAN  ANQLE  (DEGREES) 


Figure  65.  Boreelght  error  for  retloae  D. 
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APPENDIX  A 


DERIVATION  OF  EQUATION  (5.6) 


Referring  Co  Figure  21,  che  curvacure  matrix  f 4] ,  [10] of  che  incidenc 
spherical  wave  may  be  expressed  as 

a-1  0  ' 

,  •  (A.  1) 

)  a-i 

The  curvacure  matrix  of  che  inner  surface  of  che  radome  (Z.J  is 


RI1 


,-l 


(A.  2) 


From  che  knowledge  of  chese  two  curvature  matrices,  che  curvature  matrix  of 
the  refracted  ray  1-2  may  be  expressed  as  [A] 

E, 


Q1  "  +  (cos  -  cos  eJ/n)Q  ^(ej)"1 

where  0*  and  9^  are  coordinace  transformation  matrices  given  by 


(A.  3) 


cos  a. 


(A.  4) 


COS 


(A.  5) 


Simplifying  (A. 3),  we  obtain 


(A.  6) 

where  and  *r®  defined  in  (5.7). 

The  curvature  matrix  of  the  wavefront  Incident  at  point  2  la  given  by 

Q2  "  1  Wj)"1  +  b  * I*1  (A.  7) 

where  I  is  the  identity  matrix.  The  curvature  matrix  of  the  transmitted 
wavefront  at  2  can  be  calculated  similarly  to  that  at  1.  The  final  result 
is 


where  and  *22  mr*  in  (5.7). 
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